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Hadronic Light-by-Light and the Pion Polarizability 



(N 

o 

(N 
C 

Hi 

m 



43 ■ 

Oh' 
i ' 

43 



(N 
> 
ON 

o 

oo 
q 

o 

(N 



X 



Kevin T. Engcl, 1 Hiren H. Patel, 2 and Michael J. Ramsey-Musolf 3 

'California Institute of Technology, Pasadena, CA 91125 USA 
2 Department of Physics, University of Wisconsin, Madison, WI 53706, USA 
s Department of Physics, University of Wisconsin, Madison, WI 53706, USA 
and Kellogg Radiation Laboratory, California Institute of Technology, Pasadena, CA 91125 

We compute the charged pion loop contribution to the light-by-light scattering amplitude for 
off-shell photons in chiral perturbation theory through next-to-leading order (NLO). We show that 
NLO contributions are relatively more important due to a fortuitous numerical suppression of the 
leading-order (LO) terms. Consequently, one expects theoretical predictions for the hadronic light- 
by- light (HLBL) contribution to the muon anomalous magnetic moment, a^ LBL , to be sensitive 
to the choice of model for the higher momentum-dependence of the LBL amplitude. We show 
that models employed thus far for the charged pion loop contribution to aJJ LBL are not consistent 
with low-momentum behavior implied by quantum chromodynamics, having omitted potentially 
significant contributions from the pion polarizability. 



The anomalous magnetic moment of the muon, a M = 
(<7 M — 2)/2, continues to be a quantity of considerable in- 
terest in particle and nuclear physics. The present exper- 
imental value, af v = 116592089(63) x 10" 11 obtained by 
the E821 Collaboration[lH3 differs from theoretical ex- 
pectations by 3.6cr assuming the Standard Model (SM) 
of particle physics and state-of-the-art computations of 
hadronic contributions, including those obtained using 
data on er(e + e~ — > hadrons) and dispersion relation 
methods: a™ = 116591802(49) x lO" 11 (for recent re- 
views, see Ref. 0, [H[ as well as references therein). A 
deviation of this magnitude can be naturally explained 
in a number of scenarios for physics beyond the Stan- 
dard Model (BSM) , including (but not limited to) super- 
symmetry, extra dimensions, or additional neutral gauge 
bosons [6j-|8[ . A next generation experiment planned 
for Fermilab would reduce the experimental uncertainty 
by a factor of fourQ. If a corresponding reduction in 
the theoretical, SM uncertainty were achieved, the muon 
anomalous moment could provide an even more powerful 
indirect probe of BSM physics. 

The most significant pieces of the error quoted above 
for a^ M are associated with the leading order hadronic 
vacuum polarization (HVP) and the HLBL contributions: 
£o£VP( L q) = ±42 x 10 -n an d <5 a BLBL = ±26 x lO" 11 
[lfj| (other authors give somewhat different error esti- 
mates for the latter 14-23] ; but we will refer to these 
numbers as points of reference; see 23| for a review). In 
recent years, considerable scrutiny has been applied to 
the determination of a BVP (LO) from data on <j(e + e~ — > 
hadrons) and hadronic r decays. Use of the latter indicat- 
ing a somewhat smaller discrepancy between the SM and 
experimental values for a M than quoted above. Clearly, 
a significant improvement in this determination will be 
needed if the levels of theoretical and future experimental 
precision are to be commensurate. 

Here, we concentrate on the a BLBL , focusing in partic- 
ular on the contributions from charged pion loops. Sub- 
sequent to the first results from the E821 Collaboration, 



the theoretical community devoted substantial effort to 
refining the predictions for pscudoscalar "pole" contri- 
butions, which appear at leading order in the expansion 
of the number of colors Nc and which are numerically 
dominant. However, the error quoted for the charged 
pion loop contributions, which enter at subleading order 
in Nc, is now comparable to the uncertainty associated 
with the pseudoscalar pole terms. Thus, we are moti- 
vated to revisit the former as part of the effort to improve 
the level of confidence in the theoretical SM prediction 
for a BLBL . 

As a first step in that direction, we have computed the 
HLBL scattering amplitude for off-shell photons to NLO 
in Chiral Perturbation Theory (%PT). xPT is an effec- 
tive field theory for low-energy interactions of hadrons 
and photons that incorporates the approximate chiral 
symmetry of quantum chromodynamics (QCD) for light 
quarks. Long-distance hadronic effects can be computed 
order-by-order in an expansion of p/A x , where p is a typ- 
ical energy scale (such as the pion mass m n or momen- 
tum) and A x = AttF^ ~ 1 GeV is the hadronic scale with 
= 93.4 MeV being the pion decay constant. At each 
order in the expansion, presently incalculable strong in- 
teraction effects associated with energy scales of order A x 
are parameterized by a set of effective operators whose 
coefficients - "low energy constants" (LECs) - are fit to 
experimental results and then used to predict other low- 
energy observables. 

%PT has been applied with considerable success to the 
analysis of a variety of hadronic and electromagnetic pro- 
cesses (for a recent review, see e.g. 



24(), making it an 



in principle appropriate and model-independent frame- 
work for investigating hadronic contributions to a^, an- 
other low-energy observable. In the xPT analysis of the 
pseudoscalar pole contributions to a BLBL , however, one 
encounters a new LEC that cannot be determined in- 
dependent of the a M measurement itself. Consequently, 
hadronic modeling is presently unavoidable if one wishes 
to predict the anomalous moment. Nevertheless, the cal- 



culable terms in x^T can be used to test or constrain 
model input, as any credible model for the LBL ampli- 
tude must reproduce behavior in the low-energy regime 
that is dictated by QCD. Indeed, the xPT computa- 
tion of the leading In 2 term in the pion pole contribu- 
tion revealed a critical sign error in earlier numerical 
computations of the pion pole contribution [l7, 18 1. The 
sub-leading In term can be obtained from a combination 
of analytic computation [22l] and a determination of the 
relevant LEC from a determination of the ir° — > e + e~ 
branching ratio [2HJ, and it can be used to further con- 
strain the model input. 

In this spirit, we have analyzed the charged pion 
loop contribution to the LBL amplitude to NLO and 
have compared with corresponding predictions implied 
by models used in the computation of a^ LBL . The lead- 
ing order (in chiral counting) contribution is hxed entirely 
by gauge invariance and contains no unknown constants. 
As we show below, this contribution is fortuitously sup- 
pressed. As a result, higher order contributions arc likely 
to be relatively more important than one might expect on 
general grounds, rendering this quantity more susceptible 
to model-dependent uncertainties. Thus, it becomes all 
the more important that any model used for the charged 
pion contribution to a BLBL respect the requirements of 
QCD at NLO in the low-momentum regime. In this re- 
spect, we find that models utilized to date have omitted 
a potentially significant contribution associated with the 
pion polarizability, leading one to question the reliability 



of the presently-quoted value for a 



HLBL 



Below, we pro- 



vide details of the calculation leading to this conclusion. 

We compute the charged pion contributions to the 
LBL vertex function IP 1!/Q/3 through NLO from the di- 
agrams in Figure 1, expanding the result as a power se- 
ries in the external (photon) momentum and pion mass. 
The LO amplitude that corresponds to a pure scalar 
QED calculation for point-like charged pions follows from 
Fig. 1(a) and yields a finite result that is free from 
any LECs. The result contains two 0(p 4 ) structures 
that can be expressed in terms of two dimension eight 
(d = 8) operators, 32 of ] = (F 2 ) 2 = (F^F^f and 
8C^ 8) ee F 4 = F af} FPiF lX F Xa , whose coefficients arc 
given in Table U (the operators are defined to absorb sym- 
metry factors) . Naively, one would expect the magnitude 
of the coefficients to be set by 1/ (47r) 2 x 1 jm% . However, 
we find that each operator contains an additional sup- 
pression factor of 1/9 and 1/45, respectively. Thus, we 
anticipate that the NLO contributions from the graphs 
of Fig. l(b-d) will be relatively more important. 

The graphs in Figures l(b-d) correspond respectively 
to the propagator, vertex, and polarizability corrections. 
The first two classes are divergent and require the in- 
troduction of counterterms from the 0{p A ) chiral La- 
grangian. Wc carry out the calculation using dimensional 
rcgularization in d = 4 — 2e dimensions and define the 



(a) LO 7r^ contribution to LBL 



X> 1 X^" 



(b) NLO propagator correction 



(c) NLO 77T7T vertex correction 



Xx 



(d) NLO 777771 vertex correction 



FIG. 1: Representative diagrams for charged pion loop con- 
tributions to the LO and NLO to LBL amplitude. 



counterterms to remove the contributions proportional 
to 1/e — 7 + ln47r + 1 as is the standard convention 
for xPT[24j. We find that the explicit dependence on 
the counterterms needed for renormalization of the pion 
propagator is cancelled by charge and mass renormaliza- 
tion, leaving only a dependence on the 0(p 4 ) operator 
associated with the charge radius of the pion: 



£9 = 



iea 9 F„ v Tr (Q [D"E, D v tf 



]) 



(1) 



where Q = diag(2/3, — 1/3) is the electric charge ma- 
trix and E = exp(ir a n a /F n ) with a = 1, 2,3 giving the 
non-linear realization of the spontaneously broken chiral 
symmetry. After renormalization, one has for the square 
of the pion charge radius 



12 

F3 



1 

A2" 

X 



In I ^rr ] - 1 

rat 



(2) 



where the superscript "r" indicates the finite component 
after the subtraction of 1/e — 7 + In 4-/T + 1 term is per- 
formed. Choosing fi = m p and taking the experimental 
value for r 2 gives a r 9 (m p ) = (7.0 ± 0.2) x 10~ 3 for two- 
flavor xPT at C(p 4 ). Within error bars, this result is the 
same as obtained in Ref. [26| for the three-flavor case. 

The 7T7T77 vertex correction shown in Fig. 1(d) is finite, 
but the polarizability amplitude nevertheless receives an 
additional finite contribution from Cg and 



£10 = e 2 ai F 2 Tr (QEQE t ) 



(3) 
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The corresponding combination entering the LBL ampli- 
tude is «g + a r w . As the sum of the one- loop polariz- 
ability sub-graphs is finite, this combination of LECs is 
independent of the renormalization scale. An experimen- 
tal value (a r g + a r 10 ) cxp = (1.32 ± 0.14) x 1(T 3 has been 
obtained from radiative pion decay [I?} ■ As a cross check 
on the extraction of these LECs we also consider the de- 
termination of a\ from semilcptonic r-decays given in 
Ref. [28j]. Converting from three- to two-flavor %PT we 
obtain a r w (m p ) = —(5.19 ± 0.06) x 10 -3 , in reasonable 
agreement with the determination of ag(m p ) from the 
pion form factor and (ag+a\ ) from pion radiative decay. 
The resulting prediction for the pion polarizability [29j . 
which we confirm by taking the on-shcll photon limit 
of our off-shell tt + it~jj computation, disagrees with the 
latest experimental determination [3^ by a factor of two. 

The final NLO results for the LBL amplitude are sum- 
marized in Table I. To lowest order in external momenta, 
the only change from LO are polarizability corrections 
which modify the coefficient. To see the full impact 
of the (higher momentum) NLO terms, we expand our 
result to 0(p 6 ), introducing a complete basis of seven 
d = 10 four-photon operators: 
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The coefficients of these operators are given in Table II. 
At this order, both vertex and polarizability corrections 
modify the LO result. 

To obtain a sense of the numerical impact of the two- 
loop corrections, including those involving oig + a r Wl we 
utilize the values of the LECs discussed above. In the 
case of , the NLO (two-loop) contribution represents 
a ~ 20% correction to the LO term, substantially larger 
than the ~ m^/A^, ~ 0.01 magnitude one might ex- 
pect from power counting arguments. In the case of the 
d = 10 operators, the NLO corrections range from a few 
to ~ 30%. The largest impact of the charge radius cor- 
rections is on O^ 10 ' 1 (~ 30%) while the most important 

effect of the polarizability is on O^ ^ (~ 10%). As we 
discuss below, the numerical impact of the various NLO 
contributions on the low-momentum HLBL amplitude - 
while illustrating their relative importance due to the LO 
suppression - may not be indicative of their impact on 
the a^ LBL . Indeed, previous experience with the inclu- 
sion of the pion form factor in earlier work [l2 . 13. 15. 16| 



nounced than implied by these low-momentum compar- 
isons. 

TABLE I: Coefficients of lowest dimension (d = 8) operators 
contributing to the HLBL amplitude, scaled by (47r) 2 m*/e 4 . 
Second and third columns give LO and NLO contributions in 
xPT, while final column indicates the VMD result 
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TABLE II: Coefficients of d — 10 operators Oi 10 ' contributing 
to the HLBL amplitude, scaled by (47r) 2 m|/e 4 . First column 
denotes operator index n. Second and third columns give LO 
and NLO contributions in xPT, while final column indicates 
VMD result. Identifying r 2 = 6/My (see text) implies agree- 
ment between the two-loop xPT and VMD predictions for the 
charge radius contribution. 
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suggests that the effect on a BLBL may be even more pro- 



We now compare the explicit NLO results in xPT 
with the corresponding expectations for the operators 
in Tables I and II derived from models used to com- 
pute the charged pion loop contribution to a BLBL . For 
concreteness, we focus on the extended Nambu-Jona- 
Lasinio (ENJL) model adopted in Ref. [Hj]. In that 
work, the point-like contributions to the LBL vertex 
function fi"" 1 ^ arc modified by the inclusion of vec- 
tor meson dominance (VMD) type propagator functions 
V pX {k 2 ) = (g p xM 2 ~ Pp p x )/(M 2 -p 2 ) as 

^V^( Pl )V urT (p 2 )V ap ( P3 )V Pll (p,)U x ^ , (4) 

with the "vector meson mass" My in general a function 
of the photon momentum p 2 . The Ward identities im- 
ply that the p^px terms do not contribute to the overall 
LBL vertex function; hence, the replacement of Eq. (j4j is 
equivalent introducing a VMD form factor for each pho- 
ton when My is taken to be a constant. The correspond- 
ing prediction for the charge radius is (r^)vMD = 6/My. 
For My = m p , one obtains a value for r 2 in good agree- 
ment with experiment. An analogous treatment using a 
Hidden Local Symmetry approach 15, 1(| agrees with 
the ENJL prescription to 0(p 6 ). 

Expanding the right hand side of Eq. ^ to first order 
in p 2 1 M v we obtain the model prediction for the NLO 
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operator coefficients given in the last column of Table I. 
Identifying 6 /M v with the corresponding quantity that 
gives the pion charge radius , we observe that the VMD 
model reproduces some but not all of the physics that 
one expects at NLO for the LBL amplitude. In par- 
ticular, the polarizability contributions to as well 
as 5 are absent from the VMD prescription. As a 
point of principle, the results of this comparison imply 
that the VMD-type models employed for a BLBL are not 
fully consistent with the strictures of QCD for the low- 
momentum behavior of U> JjVa @ and that use of a more 
consistent model prescription is warranted. 

On a practical level, given the relative magnitudes of 
the olq + a r l0 and ctg , one has reason to suspect that the 
omission of the polarizability contribution could have nu- 
merically significant implications for a BLBL . As discussed 
earlier, a comparison of the low-momentum LO and NLO 
contributions to the low-momentum HLBL amplitude in- 
dicates that the both the charge radius and polarizability 
contributions that appear at NLO can generate substan- 
tially larger corrections than one might expect based on 
power counting, due to the fortuitous numerical suppres- 
sion of the LO terms. Moreover, the charge radius and 
polarizability contributions can have comparable magni- 
tudes in the case of some operators, while for others, one 
or the other dominates. 

At this point, one may only speculate as to the ef- 
fect on a BLBL of the previously neglected polarizability 
contribution. Nevertheless, it is instructive to refer to 
existing model computations that introduce a pion form 
factor at the 7r + 7r~7 vertices. In the original compu- 
tation of Ref. [lH, inclusion of the form factor via a 
VMD prescription reduced the magnitude of the charged 
pion loop contribution to a BLBL by a factor of three from 
the scalar QED/point-likc pion result. The subsequent 
computation using the HLS procedure yielded an even 
stronger suppression (a factor of ten) [TBI, [l6| . The ENJL 
calculation of Ref. [13[ leads to a result that is about four 
times larger than the HLS computation, but still strongly 
suppressed compared to the point-like pion/scalar QED 
limit. In all cases, the use of a VMD type procedure that 
matches onto the r\ terms for the HLBL amplitude at 
low-momentum has a much more significant numerical 
impact on a BLBL than the low-momentum comparisons 
would suggest. Given that the latter already indicate a 
substantial contribution from the pion polarizability, it 
appears important to include the corresponding physics 
in modeling the charged pion contribution to a BLBL . An 
effort to do so will be reported in forthcoming work. 
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